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We consider the stochastic dynamics near zero-temperature of the random ferromagnetic Ising 
model on a Cayley tree of branching ratio K. We apply the Boundary Real Space Renormalization 
procedure introduced in our previous work (C. Monthus and T. Garel, J. Stat. Mech. P02037 (2013)) 
. in order to derive the renormalization rule for dynamical barriers. We obtain that the probability 

7—i • distribution Pn {B) of dynamical barrier for a subtree of n generations converges for large n towards 

' some traveling- wave P„{B) ~ P*[B — nv), i.e. the width of the probability distribution remains 

04 , finite around an average-value that grows linearly with the number n of generations. We present 

numerical results for the branching ratios K — 2 and K = 3. We also compute the weak-disorder 
' expansion of the velocity v for K = 2. 

Si 

; I. INTRODUCTION 

The stochastic dynamics of classical disordered spin systems turns out to be extremely slow in the whole low- 
P5 ■ temperature phase (see for instance the books [H [1] and references therein) . The reason is that the system tends to 
' remain trapped in valleys of configurations on various scales. Within the point of view of dynamical simulations, this 
. problem has been called the 'futility' problem [Sj : the number of distinct configurations visited during the simulation 
• ^-H ' remains very small with respect to the accepted moves. The reason is that the system visits over and over again the 
^ , same configurations within a given valley before it is able to escape towards another valley. 

As a consequence, a natural idea is to formulate some appropriate renormalization procedure for dynamical barriers. 
For random walks in random media. Strong Disorder Renormalization rules have been formulated in real space 
^ ' 1^. The direct generalization of this approach to many-body systems leads to Strong Disorder Renormalization 
'"^ ] formulated in configuration space 0, because the renormalization concerns the master equation of the dynamics 
^ I defined in configuration space. These Strong Disorder RG approaches are perfect to describe correctly the hierarcal 
O ' organization of valleys within valleys. However, from a numerical point of view, since the size of the configuration 
space grows exponentially with the number of degrees of freedom of the many-body system, this approach can be 
applied numerically only for small sizes 0- 
, In a recent work [8j, we have thus introduced a different renormalization approach formulated in real space : 
K*" ■ using the standard mapping between the detailed-balance dynamics of classical Ising models and some quantum 
' Hamiltonian, we have derived appropriate real-space renormalization rule for this quantum Hamiltonian. We have 
22 solved explicitly the renormalization flow for the random ferromagnetic chain ,8], for the pure Ising model on the 
' Cayley tree ,8], and for the hierarchical Dyson Ising model [9]. In the present paper, we study the case of the random 
ferromagnetic Ising model on the Cayley tree. 

The paper is organized as follows. In section |IT1 we explain how the Real Space Renormalization approach of Q 
can be applied to the stochastic dynamics of the random ferromagnetic Ising model on the Cayley tree. In section 
I mil we derive the renormalization rules for dynamical barriers near zero temperature. In section IIVI we study the 
• • renormalization flow for the probability distribution of dynamical barriers. The cases of branching ratio K = 2 and 
. 5^ • K = 3 are discussed respectively in sections |V] and I VII with numerical results on the traveling- wave statistics of 
dynamical barriers. Section [VIII summarizes our conclusions. 



II. REAL SPACE RENORMALIZATION APPROACH 



A. Model and notations 



We consider the random ferromagnetic Ising model with the classical energy 

UiC) = -J2j^,S,S, (1) 

i<j 

defined on a Cayley tree of branching ratio K with N generations, and with free boundary conditions on all the 
boundary spins. The coordinence of non-boundary spins is thus (K + 1). The couplings Jij are independent random 
positive variables drawn with some law p(J). 
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The stochastic dynamics is defined by the master equation 

dPt (C) ^ 



dt 



Pt (C) W (C ^C)~Pt iC) Wont (C) 



(2) 



that describes the time evolution of the probabihty Pt{C) to be in configuration C at time t. The notation W (C C) 
represents the transition rate per unit time from configuration C to C, and 



(3) 



represents the total exit rate out of configuration C. We will focus here on single spin- flip dynamics satisfying detailed 
balance 



W{C^ C) = e-f^^'^^'^W {C ^ C) 



where the transition rate corresponding to the flip of a single spin Sk reads 



(4) 



(5) 



G'^'^'^'^[hk] is an arbitrary positive even function of the local field hi^ — '^^^j. JikSi. For instance, the Glauber dynamics 
corresponds to the choice 



Glauber 



1 



2cosh(/3/i) 



(6) 



B. Associated quantum Hamiltonian 

As is well-known, the master equation of Eq. [2] with the choice of Eq. [5] can be mapped via the similarity 
transformation Pt{C) = e~^'^^'^^ipt(C) onto a Schrodinger equation with the following quantum Hamiltonian involving 
Pauh matrices (ct^,ct^) @,[i3-[i3 



n 



J2 G'" ' 



(7) 



Note that in the high-temperature limit /3 — > 0, the quantum Hamiltonian of Eq. [7] for the Glauber dynamics of 
Eq. ini reduces to the standard transverse- field Ising model 



N 



n 



N - 



fc=i 



(8) 



When the coupling are random, the low-energy physics is then well described by the Strong Disorder RG procedure 
valid both in one dimension [17] and in higher dimensions d > 1 [l8l - l20j (see [2l| for a review). However here we are 
interested into the opposite limit of very low temperature where /3 — >■ -l-oo, where one cannot linearize the exponentials 
in the quantum Hamiltonians of Eq. [T) We have explained in Q how to define appropriate real-space renormalization 
rules for this type of quantum Hamiltonian in the opposite limit of very low temperature. In the following, we recall 
the main idea for the case of the Cayley tree geometry. 



C. Reminder on the Boundary Real Space Renormalization on a Cayley tree 

For a model defined on a Cayley tree, it is natural to define a boundary renormalization procedure in order to keep 
the tree topology unchanged, and to obtain recurrence equations. For instance for the Random Transverse Field Ising 
Model of Eq. El this idea has been used either within the Quantum Cavity Approach [23 - |23 | or within the Boundary 
RG approach [25( 1 . 
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FIG. 1: Basic step of the Boundary Real Space Renormalization for a Cayley tree of branching ratio K = 2: the quantum 
Hamiltonian of Eq. [9]involves 7^ + 1 = 3 spins, namely the two renormalized boundary spins (Si, &) (that represent sub-trees) 
whose renormalized dynamics is described by renormalized amplitudes (Gi, G2), and their common ancestor S, whose dynamics 
is still described by the initial amplitude G^"^[h]. The further ancestor Sa is not taken into account as a quantum spin, but 
only through the external field Ba = JaSa seen by the spin S. 



For the Hamiltonian of Eq. [71 a boundary spin Si is connected to a single ancestor spin S via some random positive 
coupling Ji > 0, so the absolute value of its local field hi = JiS takes the single value J;. As a consequence, the function 
G'^™{hi) reduces to the number G™^{Ji). As explained in [8| for the case of the pure Ising model, it is thus convenient 
to define a Boundary Real Space Renormalization as follows. The basic renormalization step concerns K renormalized 
boundary spins {Si, S2, ■■, Sk) whose renormalized dynamics is described by some renormalized amplitudes Gi (which 
are numbers and not operators) and their common ancestor spin S whose dynamics is still described by the initial 
amplitude G™*[/i] involving also the external field Ba = JaSa induced by its next ancestor spin Sa- So we have to 
study the following effective Hamiltonian for these {K + 1) spins (S'l, .., Sk: S) 



K+l 



G'' 



■ K 



Ba 



(9) 



i=i 



The physical meaning of the amplitude Gi of the renormalized boundary spin Si that represents a whole sub-tree, is 
that the largest relaxation time for this isolated sub-tree reads (see detailed explanations in Q) 



1 

2G^ 



(10) 



Near zero temperature, the largest relaxation time is simply the equilibrium time i^' needed to flip between the two 
renormalized states states Si — +1 and Si = —1 representing the two ferromagnetic ground states of the corresponding 
sub-tree, so the corresponding dynamical barrier Bi is defined by the exponential behavior 



1 

2Gi /3->-l-oo 



(11) 



After the first renormalization step, the renormalized dynamical barrier Bi are expected to grow, so that the 
renormalized amplitudes Gi will be extremely small. The most appropriate approach is then a perturbative analysis 
in the parameters Gi that may be summarized as follows (see ^ for more details) : 

(i) When = for i = 1, 2, .., K, the spins {Si, .., Sk) cannot flip and are thus frozen. So the 2^ states 



K 



Si ,82 



62 



\s> 



S=± V2cOsh/3(Etl J^S^ + Ba) 



(12) 
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are zero-energy states of Hamiltonian of Eq. |9] when Gi = for any function G™' [h] since one has 



> 



(13) 



The physical interpretation is that the spin S is at equihbrium with respect to the frozen spins (Si, .., Sk)- The other 
2^ states have a finite energy for Gi = 0. 

(ii) When the amphtudes Gi for i — 1,2, .., K are smaU, we need to diagonahze the perturbation within the subspace 
spanned by the 2^ vectors \ySt,S2,...SK ^ ^ j^^j^ f^^. 

an eigenstate via the hnear combination 



Si=±,S2 = ±,..,Sk=± 



(14) 



The eigenvalue equation — [Hk+i — A)|ua > reads 
2cosh/3(E,^,J,^,+i?,, 



1=1 

K 



2cosh/3(^f^i Jj5j+Ba) 



(S*!, , Sk) 



2cosh/3(E,^, J,^,+i3„) 



2cosh/3(E,_^, J.S*, - + S,)^2cosh/3(Ej_^, + J^S, + B,) 



=?Ai(S'i, .., -S'i, S'k) (15) 



(iii) As explained in detail in 's'l, one obtains that near zero-temperature, the Hamiltonian of Eq. [5] can be 
renormalized onto the single spin effective Hamiltonian 



for the renormalized spin (describing the full ferromagnetic states) 



(16) 



K 



\Sr - + > = 



> 



n i^j- =+>ys=+> 

\{\Sj = ->\ \s = -> 



with the renormalized amplitude (using Eq[ 



Gr = 



Ai 



2 cosh/3 Jq 



(17) 



(18) 



where Ai is the first non-vanishing eigenvalue of the system of Eq. 1151 

In the pure case studied in [s?], i.e. when all couplings Ji take the same value J, and where all amplitudes Gi take 
the same value G, we have solved Eq[TS]by taking into account the symmetry between the K branches. Here in the 
disordered case, the K branches are not equivalent anymore, but we can nevertheless derive an explicit renormalization 
rule for Gr, as explained in the following section. 



III. RENORMALIZATION RULE FOR DYNAMICAL BARRIERS 

A. Eigensystem for Ai 

To see more clearly the meaning of the eigensystem of Eq. 1151 it is convenient to introduce 



Tq{Si, .., Sk) = 



\ 



K 



2cosh/3(^ J,5, + B,) 



which solve the system of Eq. [15] for A = 0, and the ratios 



A{Si,..,Sk) 



T\i{Si, ..,Sk) 
To{Si, .., Sk) 



(19) 



(20) 
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so that Eg [TSl reads 







K 



Eg 



Tq (5*1, ., Si, ., Sk) 



A 



K 



T§{Si,., Sj — 0, ., Sk) 
ToiSi, ., 5^, ., Sk) 



A(5i,..,-^„..,5k)(21) 



To compute the lowest non-vanishing eigenvalue Ai of Eq. [151 it is consistent to set Ai = in all equations except 
at the two extreme cases where all spins have the same value, either Sa or —Sa, where Ba = JaSa- The two extreme 
components should be orthogonal to Eq and thus read 



Tx,{Sa,Sa,...,Sa) = To{- Sa, - Sa, .. - Sa) = . 2 COsh J, - J,) ~ 6 ^ (^^-^ -'^") 



( Sa, Sa, Sa) — To[Sa, Sa, Sa) — 

i.e. the two extreme values of the ratios of Eq. [50] 



K 



2cosh/3(^ J, + Ja) 



(22) 



■A{Saj Sa, Sa) 



To{~Sa, ^Sa, -Sa) \/ ^ ^^'^ /3(E 1 1 J^ " Ja 



A[ Sa, Sa,---, Sa) — 



To{Sa, Sa, Sa) 

Toi^Sa, Sa, Sa) 
To{ Sa, Sa, ~Sa) 



2 cosh I3ij:f^^ J, + Ja 



2cosh(3{j:f^, J, + Ja) 
2cosh(3{j:f^,J,-Ja) 



(23) 



satisfy Eq[5T]with Ai. So Ai can be computed near zero temperature as 



Ai =X^G, 

i=l 
K 

-Eg. 



T^Q^Sa, Sa, Sa, Si — 0, Sa, Sa, Sa) A^Sa, Sa, Sa, Si — Sa, Sa, Sa, Sa) 



T^^Sa, Sa, - Sa) 



1 - 



A(^Sa, Sa, Sa) 



A{Sa, Sa, - Sa, Si — Sa, Sa, Sa, - Sa) 
A{Sa, Sa, Sa) 



(24) 



whereas all non-extreme values satisfy Eq[5T]with Ai = 0, so we may drop the common denominator T^{Si, ., ., Sk) 
to obtain 







K 



E Gj^o •> Si — 0, ., Sk) 



,i=l 



K 



A[Si, , Sk) — E G^Tq (5i, ., Si — 0, ., Sk)A{Si, .., —Si, .., Sk) 



(25) 



In the disordered case, we expect that the dynamical transition between the two ferromagnetic states will be dominated 
by a single dynamical path near zero temperature. So let us now compute the dynamical barrier associated to a given 
dynamical path. 



B. Dynamical barrier associated to a given dynamical path (1,2, --,K) 
In this section, we consider the given dynamical path called (1, 2, .., K) 

A(y Sa, —Sa, Sa, Sa, —Sa) ^ A{^Sa, Sa, —Sa, —Sa, Sa) ^ A(^Sa, Sa, Sa, Sa, —Sa) 

... O A(^Sa, Sa, Sa, —Sa) ^ ^(5*0, 5*0, Sa, Sa, Sa) (26) 

Then we are left with a one-dimensional problem with the notations 

(7(0) = A{-Sa, -Sa, -Sa-, Sa) ^ -e^'^ 
q[l) ^ A{Sa,-Sa,-Sa,—,-Sa) 
q{2) = A[Sa, Sa, —Sa, —Sa) 

q{K-l) =A{Sa,Sa,-,Sa,-Sa) 

q{K) ^AiSa,Sa,-Sa,Sa)^e-^-'^ (27) 



Eg [24l reduces to 

and Eq. becomes ior 1 < k < K — 1 

q{k)=p^{k)q{k-l)+p+{k)q{k + l) 

with the probabilities 

GkTQ {Sa, Sa, Sk — 0, —Sa, ^Sa, ~Sa) 



P-{k) 
P+ik) 



GkTQ^Sa, Sa, Sk — 0, —Sa, —Sa) + Gfc+iTg (^a, , 5*0, S'/c+i — 0, — Sa, —Sa) 

Gk+lT^i^Sa, Sa, Sk+1 — 0, —Sa, , Sa) 
GkT^i^Sa, , Sa, Sk — 0, Sa, Sa) + Tq (6*0 , , Sa, Sk+1 = 0, Sa, Sa) 



normalized to p^{k) +p^{k) = 1. 

The solution can be obtained by [1, [26 



^^^ Ro{k,K-l) , RK{l,k) 
q{k) = g(0) „ ,^ -T + q{K)- 



RoiO,K-l) 'Rk{1,K) 



using Kesten variables [27 



Ro{K,K-l) =0 
i?o(-?f - 1) =1 



'-^ p~(n) 



and 



Rk{1,0) -0 
i?A-(l,l) =1 

k—1 m / \ 

i?.(i,fc>2) 

7n— 1 n— 1 -1 I ^ ^ 
K-1 



p n K) -1 I V ttp-(") -1 I P-(^) I , p-(i)p-(2)...p_(jr-i) 

' ^ i'lMp+W %+(l) %+(l)p+(2)^-^P+(l)P+(2)...p+(i^-l) 

Plugging this solution into Eq [28] yields 



g(i^)i?o(0,i^-l) 

1 + e^^''^ 



RoiO,K-l) 
so the renormalized amplitude of Eq. [TS] becomes 



2 cosh ^Ja i?o(0,-f4r- 1) 

So using Eq. [T2| we need to compute the ratios of the two probabilities of Eq. [50] 

P+{k) ^ Gk+lT§{Sa, Sa, Sk+1 = 0, -Sa, Sa) ^ Gfc+i 2 COsh 1 J^ ^ Etfc+2 J^ + Ja) 
p.{k) GkT^iSa, Sa,Sk = 0, -Sa, Sa) Gfc2 COsh /3(Etl' J^ " Etk+1 + Ja) 
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as well as the products 



n 



GKTQ{Sa, Sa, Sa + , Sk — 0) 



P+jk) ^ 

p^{k) G,nTiiSa,Sa, Sa, S,n = 0, Sa, ~Sa, Sa) G„2 COSh /^(E:^^^ J, - E.=™+1 + Ja) 

Our conclusion is thus that the dynamical path (1,2, ..,i4r) has a renormalized amplitude Q^^'^'- '^^ given by 



(37) 



1 



G 



(1,2,. 
R 



Gk 



G 



K 



' Gj^2cosh/3(^t7^ J, + J,) 

G,„2cOsh/3(X;-lT^ Ji ~ Ejirn+l + Ja) 



1 + E 



g/3(JK-Jo)2cosh/3(^ J, + J„) ^ 



\ G,„2cosh/3(X;r=/ - E^m+l J^ + 



if 

E 



-1 ^m<^ 

In terms of the dynamical barriers near zero temperature introduced in Eq. [TT] 



Gi oc e 

,3-> + oo 



(38) 



(39) 



Eq. [551 yields that the renormalized barrier 
ature is given by 



(1,2,. 



associated to the dynamical path (1,2,.., K) near zero temper- 



In 1 

(12 if) g(1.2,.../f) 

Sn.' = lim ^ = max 

" ^->-|-oo p l<m<K 



K 



max 

l<m<if 



1=1 

/" m — l 



K 



^ ^ ^ i ^ ^ Ji ~^ Ja 

K 



Bjn + Jtu - Ja + 2 min ^ Ji + Ja; ^ Jj 



z=7n+l 



(40) 



C. Optimization over the K\ dynamical paths 

We now have to consider the K\ possible dynamical paths : for a given permutation tt of the K renormalized spins, 
the dynamical barrier associated to the path (7r(l), 7r(2), .., iriK)) reads by adapting Eg l40l 



^(.(l),.(2),..,.(if)) ^ 

" l<m<K 



K 



fi7r(m) + J7r(m) ^ Ja + 2 min I ^ J.,r(i) + Ja', ^ J,r(i) 

i—m-\-l / 



, i=l 



(41) 



We now have to choose the dynamical path, i.e. the permutation tt leading to the smallest barrier. So the final 
renormalization rule is that the renormalized barrier Br is given by the minimum of Eq. 21] over the K\ possible 
permutations 



Br — mm I Op 

7r V 

/ r /m-l K 

^7r(m) + J7r(m) - Ja + 2 miu ^ Jjr(i) + Jq; ^ Jt,(^ 



mm max 

TT \ l<m<K 



. i=l 



2=m+l 



(42) 



D. Limit of the pure case 

In the pure case, all couplings Ji have the same value Jq and all renormalized boundary spins Si have also the same 
barrier Bi = B, so that the renormalization rule of Eq. 01] reduces to 

Br — B + 2 Jo max [min (to; K — to)] (43) 

l<jn<K 
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For even K, the maximum is reached for m = K/2, whereas for odd K ^ the maximum is reached for m = {K ± l)/2, 
so that Eq. H51 yields 

Br =B + KJq for even K 

^ B + (K ^ l)Jo for odd K (44) 

This recurrence thus leads to the following linear growth 

Bn — Bn=o = tiKJq for even K 

= n{K - 1) Jo for odd K (45) 

with the number n of generations, in agreement with the more detailed results of [8,] containing explicit combinatorial 
prefactors coming from the degeneracy of the K\ dynamical paths. The linear growth with the number n of generations 
is in agreement with previous works of physicists [28l - l30| and of mathematicians [sil - fs^ . It turns out that the slope 
{K — l)Jo for odd K of Eq. H51 coincides with the slope obtained in [2^, where a so-called 'disjoint strategy' is 
optimal, whereas the slope KJq for even K of Eq. H5] differs from the slope J{K — 1) obtained in "2^ [2^, where a 
so-called 'non-disjoint strategy' is optimal. We refer to Refs t28i-.32i] for more explanations on the differences between 
disjoint/non-disjoint strategies. Here it is clear that the renormalization procedure making coherent clusters of spins 
within sub-trees corresponds to the disjoint strategy. Within the disjoint strategy, the renormalization of Eq. H^for 
dynamical barriers is close to the recursions written in Refs [28r,-32j , even if there exists a difference : in our approach 
based on the quantum perturbation theory of section lll C[ the spin S is considered to be at equilibrium with respect 
to its neighbors, whereas in previous approaches [28|-t32|. the question of the time where the roots flips with respect 
to the flips of the sub-trees has been taken into account differently (see the notions of 'directed barriers', of 'anchored 
barriers', etc.). Let us also mention that the case of 'disorder in the degrees of the nodes' has been studied by Henley 
[28j . via the model of critical percolation on the regular Cayley tree. 



IV. RG FLOW FOR THE PROBABILITY DISTRIBUTION OF DYNAMICAL BARRIERS 



Recurrence for the joint distribution of dynamical barriers and couplings 



The renormalization rule of Eq. S^lfor dynamical barriers yields that the joint probability distribution Vn{B, J) of 
the barrier B and the coupling J evolves according to the following iteration in terms of the probability distribution 
p{J) of the random ferromagnetic couplings 



Pn+l(.BH, Ja) = p{Ja) \J{ j dBMPn(.B,, J,) 



Br + Ja — rnin max 

TT \ l<m<K 



K 



B 



>/7r(m) + 2 min ^ J^(i) + Ja; ^ J^ 



(0 



4=1 



i—m+1 



with the partial normalization 



j dBPn{B,J)= p{J) 



(46) 



(47) 



The initial condition for the Glauber dynamics of Eq. [6] near zero temperature G'^qI^^^^^ [h] — 2cosh(i3h) ' 
corresponds for an initial boundary spin Si of the tree submitted to the local field h — JiS to the dynamical barrier 
Bi = Ji , so that the joint distribution at generation zero reads 

Vitt'-iB,.]) ^p{J)S{B-J) (48) 

The form of Eq. US] actually shows that it is more convenient to replace the variable B hy B = B + J , so that the 
joint distribution Pn{B,J) evolves according to 



P„+i(B, Ja) = p{Ja) / dB,dJ,P„{B,, J,)^ 



B — min max 

IT \ l<m<K 



K 



B 



IT [in) 



2min ^ J^(i) + Jq; ^ J^ 



i=l 



(49) 



with the same partial normalization as in Eq. |47] 

J dBP^{B,J)^ p{J) 

and the initial condition (Eq. |48|) 

palauber ^ ^ p{J)5{B - 2 J) 



(50) 



(51) 



B. Recurrence for integrated probability distributions 



To factorize the minimum and maximum functions involved in the delta function of Eq. 1491 it is convenient to 
introduce the two complementary integrated distributions 



F^{B,J) = / dB'PniB\j) 
Jb 

Gn{B,J) = f dB'Pn{B',J)=^p{J)-Fn{B,J) 



as well as the notation 



K 



- 2 min I ^ J,(,) + Ja; ^ J^(,) 

\ i—1 2— m+1 



Then the recurrence of Eq. [15] yields 

r' + OO 



Fn+l{B,Ja) = / dB'Pn+l{B',Ja) 
JB 

= P{Ja) (n / dBMPn{B„.,.h)\ e 



i=l 
K 



B < min ( max 

TT \l<m<K 



(52) 



(53) 



p{ja) n / dBMPniB,, Ji)\ ( ^ ^ ^^^^K r"(") + ^- 



T(m) 



(54) 



Let us consider that we have ordered the K\ permutations tt, so that the last product of Eq. [S3]can be expanded 



as 



B < max 

l<m<K 



n 1 



B > max 

l<m<A' 



(m) 



Ki. 



B > max 

l<m<K 



-Oir,(m) -I- '^7r,(m) 



i+E(-ir E n 



(55) 



P=l {7ri<7r2..<7rp} g=l "1=1 

For each permutation tt^, we may replace the product over m = 1, 2, .., K, by the product over k = iTq{m) to obtain 

T(m) + "^^^.^ 



B < max 

l<m</i' 



P=l {iri<7r2..<7rp} g=l fe=l 

if! if p 

i+E(-i)'' E nn^[ 

P=l {iri<7r2..<7rp} /c=l <3=1 

if! if 

i+E(-i)'' E n 

P=l {7ri<7r2<..<7rp} /£=1 



B>Bk + 



B>Bk+ max S'f'''^ 
i<«<p 



(56) 
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Then Eg [54l becomes 

K 



F„+i (B, Ja) = p{Ja) { n / dBd-hPniB.,, J, 



*=!■' / |_ P=l {7ri<7r2..<7rp} fc=l 



p=l {Tri<7r2<..<irp} \fc=l ' 

K\ / K 



B> Bk+ max ^''"^ 

~ 1<9<P 



Bk<B~ max S* , 
i<«<p 



P=l {7ri<7r2<..<Trp} \fc=l ^ " / 



So we obtain the following closed recurrence for the integrated probability Gn of Eq. [52] 

p=l {7ri<7r2..<Trp} 1<9<P y 

with the initial condition of Eq. [?T] 

G'^^f^^iB,J) ^ dbP^'^^''^^ib,J) = piJ)9iB-2J) =0 for B<2J 

— OO 

= p( J) for B>2J (59) 

C. Traveling-wave solution 

For disordered models defined on Cay ley trees , it is very common to find that the probability distribution P„(A) 
of some observable A propagates with a speed v and a fixed shape P* as the number n of generations grows 

P„(A) ~ (60) 

This means that the average value grows linearly with n, whereas the width around this averaged value remains 
finite. This property was discovered by Derrida and Spohn (ssj on the specific example of the directed polymer in a 
random medium, where the observable A of interest is the free-energy, and was then found in various other statistical 
physics models defined on Cayley trees [36} . This traveling- wave propagation of probability distributions have also 
been found in quantum models defined on Cayley trees, in particular in the Anderson localization problem (37l - l4]| . 
The conclusion is thus that the recursion relations that can be written for observables of disordered models defined 
on trees naturally lead to the traveling wave propagation of the corresponding probability distributions. 

So here, it is natural to expect that the solution of Eq. [51] starting from the initial condition of Eq. [55] will be a 
traveling-wave with some velocity v 

Gn{B,J)~G*{b = B-nv,J) (61) 
where the stable shape G* {b, J) of the front satisfies the equation 

G*{b - V, Ja) = -p{Ja) E(-l)'' E n / - S'^u'\jk) (62) 

and the following boundary conditions at infinity 

G*(6,J) ~ 

G*(&,J) ^ p{J) (63) 

This means that a barrier i?„ for a sub-tree of n generations reads 

B„ = nv + b (64) 

where & is a random variable of order 0(1). In the following sections, we have checked numerically for ivT = 2 and 
K = 3 that the statistics of dynamical barriers indeed follow the traveling-wave form of Eq. [HI] 
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V. CASE OF BRANCHING NUMBER K = 2 
A. Renormalization rule for dynamical barriers 

Here there are only K\ ~ 2 permutations tti = (1,2) and 112 — (2,1), so the barriers associated to these two 
permutations are (EqHT|) 

Sjj^'^^ = max [Bi + Ji-Ja + 2 min( J2, Ja); B2 + J2- Ja] 

S^''^^ =max[62 + J2- Ja + 2 min(Ji, J,); Si + J1-JJ (65) 
and the optimization of Eq. |42] reads 

Br + Ja = min (max [Bi + Ji + 2 min( J2, Ja); B2 + J2] ; max [B2 + J2 + 2 min( Ji, Ja)]Bi + Ji]) (66) 

B. Recurrence for the integrated probability 

For K = 2,we have to consider the two permutations tti — (1, 2) and 112 — (2, 1), so the variables of Eq. [551 reads 

5^"^^ -2min(J2,J,) 
4"'' = 
^^"^^ = 

4"'' =2min(Ji,Ja) (67) 

and Eq [58] becomes 

G„+i(B, Ja) = p(Ja) y dJi y dJ2G„(S - Ji)G„(B - S^^'\j2) 

+p{Ja) I dJi J dJ2Gn{B - S'f'\ji)Gn{B ~ S^^^\J2) 

~p{Ja) J d.hj dJ2G„(i?-max(5("^\5i"^)),Ji)G„(B-max(4"^\5f^)),J2) 
= piJa) J dJi J dJ2G„(B-2min(J2, Ja), Ji)G„(B, J2) 
+p{Ja) J d.h J dJ2G„(B, Ji)G„(B-2min(Ji, Ja), J2) 

-p{Ja) J dJi j dJ2Gn{B - 2min(J2, Ja), Ji)G„(B - 2min(Ji, Ja), J2) (68) 
i.e. finaUy, the integrated probabihty G„ satisfies the recurrence 

G„+i(B, Ja) = p(Ja) J dJi j dJ2Gn{B - 2min(J2, Ja)), Ji) [2G„(B, J2) - Gn{B - 2min(Ji, Ja), J2)] (69) 

C. Traveling- wave Ansatz 

The travehng-wave Ansatz of Eq. [61] yields the foUowing equation for the stable shape G* (6, J) of the front 

G*{h-v,Ja) ^p{Ja) j d.h j dJ2G*(6-2min(J2,Ja)),Ji)[2G*(6,J2)-G*(6-2min(Ji,Ja),J2)] (70) 
with the boundary conditions of Eq. [63] 
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D. Weak-disorder expansion 



For the pure case p{J) = 6{J — Jq) discussed in section ITlI Dl the solution is simply 

G;ureib,J) ^ SiJ ~ Jo)e{b) = p{j)e{b) 

with the velocity of Eqs [44] and [45l 



'-'pure 



2Jo 



(71) 



(72) 



Let us now consider the case of weak disorder, where the disorder distribution p{J) displays a small width around 
Jq. Then we expect that the width of the barrier distribution will also be small with respect to (2Jo). In the region 
where G*(b — v, Ja) is small, the appropriate linearization of Eq. [TOlconsists in replacing G*(6, J2) by its asymptotic 
behavior p{J2) fEg [63)) and by neglecting G*{b — 2min(Ji, Jq), J2) that would give a quadratic contribution, so that 
we obtain in this weak-disorder regime 



G*{b - V, Ja) ~ 2p{Ja) J dJi J dJ2p{J2)G*{b - 2min(J2, Ja), Ji) 

This means that the integral 



g*{b) = J dJG*{b,J) 

satisfies the closed linear equation 

9*{b-v) ~2 j dJaPiJa) J dJ2p{J2)9*{b-2mm{J2,Ja)) 



The exponential shape of coefficient p 



g*{b) c e 



fib 



is then a solution of the linearization of Eq. [75] if the velocity v(p) satisfies 

e-'^''^^) = 2 J dJap{Ja) J d.hp{J2)e^^^'^'''^''^''^ 

The velocity as a function of p 



v{p) 



presents the two limiting behaviors 



1 



■In 



2 / dJapiJa) / dJ2p(J2)e-''^'""(-'='^") 



, , ln2 

^-00 

M^O p 



v{p) 



^ 2J +i^^2J 



(73) 



(74) 



(75) 



(76) 



(77) 



(78) 



(79) 



where Jmin is the minimal value of the disorder distribution p(J). Between these two limits, there exists a maximum 
value V* at p* such that v'{p) — 0, and this is this velocity v* that will be dynamically selected [12, It can be 
computed for various disorder distributions p(J). Let us now consider two explicit cases. 



1. Exponential distribution for the disorder 



For the exponential distribution 



p{J)^9{J-Jo)^e-^ (80) 
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Eq [78] reads 



.(M) =2Jo + ilnii^ 



We are looking for the point /i* where this function is maximum 

=v'{fi,) 



1 , fl + fi,A 



^l% \ 2 J ^*(i + /i*A) 

i.e. = x/A where x ~ 3.31107 is the root of the numerical equation 

'1 + a;' 



In 



l + x 



The corresponding selected velocity reads 



v^^v[pi^) = 2Jo H In I 

A, fl + x 
= 2Jo + -ln^ 



l + x 



In conclusion, we obtain that the velocity grows linearly for weak disorder A 

V* = 2 Jo + 0.23195A 



(81) 



(82) 



(83) 



(84) 



(85) 



For the box distribution of width A 



Eg [751 becomes 



2. Box distribution for the disorder 



p{j) = -e{jo < J < Jo + A) 



v{ii) = 2 Jq In 

/i 



-^A-A + 2/lA - 1 



^2A2 



The corresponding selected velocity v* — v(fi*) (where v'{fj,*) = 0) grows again linearly for weak disorder A 

V* = 2 Jo + Aw* ~ 2 + 0.206826A 



(86) 



(87) 



(88) 



E. Numerical results obtained via the pool method 

If one wishes to study numerically real trees, one is limited to rather small number n of generations, because the 
number of sites grows exponentially in n. From the point of view of convergence towards stable probability distributions 
via recursion relations, it is thus better to use the so-called 'pool method' that allows to study much larger number 
of g enerations. The pool method has been used for disordered models on trees [l^jEO, SH or on hierarchical lattices 
[44| - [4a |. The idea of the pool method is the following : at each generation, one keeps the same number Mpooi of 
random variables to represent probability distributions. Within our present framework, the probability distribution 
Pn{B, J) at generation n will be represented by a pool of Mpooi — 10^ couples {Bi, Ji). To construct a new couple 
[Bii, Ja) of generation (n + 1), one draws K couples (S^, Ji) within the pool of generation n and apply the rule of Eq. 



The numerical results obtained via the pool method for the box distribution of the disorder of Eq. [86] with Jo = 1 
are shown on Fig. [21 On Fig. [H(a), we display the convergence towards the traveling- wave form for the finite disorder 
width A = 1. On Fig. [5] (b), we show that the velocity f (A) as a function of the disorder strength has for tangent 
the weak-disorder expansion of Eq. [55] near A — )■ 0. 
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FIG. 2: Cayley tree of branching ratio K = 2 with the box distribution of disorder of Eq. [86]with Jo = 1 (a) For A = 1 : the 
probabihty distributions Pn{B) of the dynamical barrier B at generation n = 5, 10, 15, 20, 25, 30, 35, 40, 45, 50 correspond to a 
travehng-wave with a fixed shape Pn{B) P*(b = B ~ nv) (after an initial transient), (b) Velocity v of the traveling- wave as 
a function of the disorder width A and agreement with the weak-disorder expansion at first order in A of Eq. 1881 



VI. CASE OF BRANCHING NUMBER if = 3 
A. Renormalization rule for dynamical barriers 

Here there are K\ = Q permutations so the barriers associated to these permutations are fEg I4ip 





(1,2,3) 


= max 




h Ji 


-Ja^ 


- 2 min( J2 - 


H J3; Ja) 


B2 


+ J2 


-Ja 


+ 2 min( Jl - 


h Ja, J3) 


Bs^ 


h J3 


-Ja] 


BT 


(2,1,3) 


= max 


[B2^ 


h J2 


-Ja- 


- 2 min( Ji - 


1- J3; Ja) 


Bi 


+ Jl 


-Ja 


+ 2 min( J2 - 


h Ja, J3) 


Bs^ 


h J3 


-Ja] 


BT 


(1,3,2) 


= max 


[Bi^ 


h Ji 


-Ja- 


- 2 min( J2 " 


1- J3; Ja) 


B3 


+ J3 


-Ja 


+ 2 min( Jl - 


h Ja, J2) 


B2^ 


h J2 


-Ja] 


BT 


(3,1,2) 


— max 


[Bs- 


h J3 


-Ja- 


- 2 min( J2 " 


1- Ji; Ja) 


Bi 


+ Jl 


-Ja 


+ 2 min( J3 - 


h Ja, J2) 


B2^ 


h J2 


-Ja] 


BT 


(2,3,1) 


— max 


[B2- 


h J2 


-Ja- 


- 2 min( Ji - 


1- J3; Ja) 


B3 


+ J3 


-Ja 


+ 2 min( J2 - 


h Ja, Jl) 


Bi^ 


h Jl 


-Ja] 


BT 


(3,2,1) 


= max 




h J3 


-Ja- 


- 2 min( J2 - 


1- Ji; Ja) 


B2 


+ J2 


-Ja 


-|- 2 min( J3 - 


h Ja, Jl) 


Bi^ 


h Jl 


- Ja] (89) 



and the optimization of Eq. |42] reads 

Bn ^ min 4^'^'^)) (90) 

The corresponding iteration of Eq[58]for the joint probability distribution is rather lengthy and not very illuminating, 
so we vifill not write it here, but present instead our numerical results. 

B. Numerical results obtained via the pool method 

The numerical results obtained via the pool method for the box distribution of the disorder of Eq. [86] with Jo = 1 
are shown on Fig. [3] We find again that the statistics of dynamical barriers is a traveling wave (see Fig. [3] (a)). The 
corresponding velocity v{A) as a function of the disorder strength A is shown on Fig. [3] (b). 

VII. CONCLUSION 

To study the stochastic dynamics near zero-temperature of the random ferromagnetic Ising model on a Cayley 
tree of branching ratio K, we have applied the Boundary Real Space Renormalization procedure introduced in our 
previous work ^8] to derive the renormalization rule for dynamical barriers. The main outcome is that the probability 
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FIG. 3: Cayley tree of branching ratio K = 3 with the box distribution of disorder of Eq. [86]with Jo = 1 (a) For A = 1 : the 
probabihty distributions Pn{B) of the dynamical barrier B at generation n = 5, 10, 15, 20, 25, 30, 35, 40, 45, 50 correspond to a 
travehng-wave with a fixed shape Pn{B) ~ P*(6 = B — nv) . (b) Velocity v of the traveling- wave as a function of the disorder 
width A 



distribution Pn{B) of dynamical barrier for a subtree of n generations converges for large n towards some traveling- 
wave Pn{B) P*{B ^ nv), i.e. the width of the probability distribution remains finite around an average- value that 
grows linearly with the number n of generations. We have presented numerical results for the branching ratios K ^ 2 
and K = and we have computed the weak-disorder expansion of the velocity v for K — 2. 

As explained in section IIV C[ the recursion relations that can be written for observables of disordered models 
defined on trees naturally lead to the traveling wave propagation of probability distributions. So we expect that for 
other disordered statistical models defined on trees, the statistics of dynamical barriers should be also described by 
travelling- waves. 
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